Surface scattering and band gaps in rough waveguides and nanowires 
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The boundaries of waveguides and nanowires have drastic influence on their coherent scatter- 
ing properties. Designing the boundary profile is thus a promising approach for transmission and 
band-gap engineering with many applications. By performing an experimental study of microwave 
transmission through rough waveguides we demonstrate that a recently proposed surface scattering 
theory can be employed to predict the measured transmission properties from the boundary profiles 
and vice versa. A new key ingredient of this theory is a scattering mechanism which depends on the 
squared gradient of the surface profiles. We demonstrate the non-trivial effects of this scattering 
mechanism by detailed mode-resolved microwave measurements and numerical simulations. 
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The coherent scattering through systems with sur- 
face roughness is a ubiquitous phenomenon w hich oc- 
curs on vastly different length and time scales.'^ The 
effects induced by surface scattering often are the 
key for the understanding both of natural phenom- 
ena, like the scattering of underwater waves at a 
rough ocean seabed,'^ as well as of man-made devices 
like optical fibers and waveguides ,1^^ photoni c cry stal 
devices ,IHHini met amat erials thin metallic films ,'^^'^2 lay- 
ered structuresj ^^ * ^^ ! nanowires,'^^'^ in optical diffrac- 
tion tomographs,!^ and confined quantum systems in 
general.-^ For some devices, like quantum cascade 
laser^^ and gated graphene nano-ribbon^^SHSU i]^q scat- 
tering at rough boundaries was identified as one of the 
dominant factors that limits the device performance. 
Surface roughness effects might hold the key for the 
explanation of anomalously large persistent currents in 
metallic ring^^H a,nd are actively used to c ontrol gravi- 
tationally bound quantum states of neutronj ^^ l ^^ l as well 
as to enhance the thermo-electric performance of nano- 
wires.'^Sl The understanding of all these phenomena rests 
on a predictive surface scattering theory that relates the 
properties of a rough surface to the transmission charac- 
teristics of the corresponding device and vice versa. 

It is here that our analysis sets in to show that conven- 
tional surface scattering theories for the coherent trans- 
mission through waveguides lack significant ingredients. 
If the aim is to determine the transmission properties of 
a system from a given rough boundary profile, numerical 
calculations can sometimes help to ov ercom e these defi- 
ciencies of existing analytical methods .^^^l^II jf^ however, 
the aim is to solve the inverse problem, i.e., to design a 
boundary such as to obtain desired transmission proper- 
ties, improved analytical models are indispensable which 
provide the functional dependence of the transmission 



properties on the boundary profile. It is this important 
missing link that has recently been provided in terms of 
an analytical surface scattering theory (SST),'^^t2ll-^y}jjch 
we test here for the first time. 

This new SST predicts how mode-specific scattering 
lengths in waveguides depend on the details of a sys- 
tem's surface roughness. It introduces the previously 
neglected square- gradient scattering (SGS) mechanism 
and predicts that this new scattering mechanism has to 
be considered together with the conventional amplitude 
scattering (AS) mechanism.l^l Although SGS is related 
to higher-order terms in the disorder strength a it can 
even be the major scattering mechanism in systems with 
modest disorder. SST focuses, however, on long-range 
correlations, which seems to restrict its applicability to 
very long waveguides. In this Rapid Communication we 
resolve this impasse by demonstrating that the validity 
of SST extends to short, individual waveguides as well. 
We thereby provide the key for an accurate description 
of transmission through surface-disordered waveguides in 
general and a solid basis for a deeper understanding of 
many of the applications mentioned above. In addition, 
our results allow us to impose desired transmission prop- 
erties on a waveguide through the roughness of its bound- 
aries. 

To demonstrate this explicitly we perform a proof-of- 
principle experiment using the microwave setup shown in 
Fig. [1] The studied scattering system consists of a pla- 
nar waveguide with metallic walls with a "rough" central 
part of length L = 80 cm which features two undulating, 
symmetric boundaries, z ~ ±[d/2 + cr^(a;)]. The aver- 
age width of the scattering part is d = 10 cm and the 
roughness variance is a^. In the experiment we choose a 
modest disorder of a/d « 9% (ct « 8.7mm). The bound- 
ary profile function ^{x) has zero mean, (^(a;)) = 0, unit 
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FIG. 1. (Color online) Experimental setup with exchange- 
able rough boundaries. The photo shows the considered 
waveguide with the surface disordered region in the center 
and the lifted top plate. Step motors are attached to the top 
plate to shift the input and output antenna in z direction. 
Absorbers at both ends suppress reflection (indicated in the 
sketch below). 



variance, (C^(a;)) ~ 1, and features a binary correlator 
defined by 



(1) 



In the convention of the SST the angular brackets stand 
for an ensemble average over different realizations of ^(x). 
However, by assuming the ergodicity of ^(a;), the average 
can be performed along a single surface profile. There- 
fore, the quantity Wix) can also be treated as the auto- 
correlation function of the boundary profile used in the 
experiment. 

For vanishing disorder the wave vector is quantized 
along the transverse direction (z), yielding n propagating 
modes whose wave numbers in guiding direction (x) at 
frequency u are given by 




(2) 



The chosen up-down symmetry of the disorder profile 
prevents mixing of odd and even modes during the scat- 
tering process. We also restrict our study to the fre- 
quency range in which we can assume single-mode scat- 
tering for the odd and even modes. 

The localization length lI'"'^'' is defined via the fa- 
mous expression for the average logarithm of the trans- 
mittance, (lnT„) = —2L/Ln°'^\ showing an exponential 
decrease with the sample length L. In the present exper- 
iment we deal with a single, finite boundary realization. 
Thus the partial transmittances r„ correspond to the av- 
erage over the most probable (representative) realizations 
of the surface disorder ,^211 



(T„) = exp (-2i/i('-)) 



(3) 



For isolated modes the localization length Ln°'^'' of the 
nth propagating mode (n = 1,2) is twice the corre- 
sponding backscattering length, Ln°^^ = 2Ln \ a quan- 
tity which has been calculated with SST,'22l 
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Here, the two different terms originate from the ampli- 
tude scattering (AS) and the square-gradient scattering 
(SGS) mechanisms, respectively: 
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W{2kn) , (5) 
5(2fc„) . (6) 
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Correspondingly, these transmission characteristics are 
specified by two different roughness power spectra,'22l 



S{kx 



yV(a;) exp(— ifc^x)da; , (7) 
yV"2(a;)exp(-ifc^x)dx. (8) 



The first one, W{kx), is the Fourier transform of the 
roughness-height correlator ([I]), and is therefore called 
roughness-height power spectrum. The second term, 
S{kx), is referred to as the roughness- square- gradient 
power spectrum since for gaussian random processes it 
is the Fourier transform of the roughness-square-gra dient 
correlator, {C^ (x) {x')) - {^^ {x)f = 2W"^ix-x')^^ 
In order to observe the two scattering mechanisms sep- 
arately, we generate a rough boundary profile ^(x) with 
a power spectrum W{2kn) being non-zero only within 
a small fc„-range. Then the localization length and po- 
tential transmission gaps outside this range will be deter- 
mined entirely by S'(2fc„). For this purpose we choose the 
roughness-height power spectrum in the following rectan- 
gular form: 



Wik, 



2(fc+-fe_) ■ 

0, 



for 2fc„ < 
otherwise. 



kx\ < 2k+ 



(9) 



respecting the normalization W(0) — 1. In our ex- 
periment we specify the correlation parameters k^ as 
follows: /c_ = 70 m~^ and fc+ = 80 m~^. Trans- 
mission properties of rectan gular spectra ^ have 
been investigated thcoreticalljP^Ell gj^d tested in mi- 
crowave experiments'^^ '^^ in one-dimensional and quasi- 
one-dimensional disordered systems with discrete scat- 
terers. 

The disorder ^(x) with a desired power spectrum 
W{kx) was constructed with a widely used convolution 
method (see, e.g., Ref . IMl and references therein). First, 
starting from a prescribed power spectrum Wlk^), we de- 
rive the modulation function G{x) whose Fourier trans- 
form is yJW{kx), 



G{x) = 



2-K 



Vmk ) exp {ikxx) . (10) 



In the second step the disordered profile ^{x) is generated 
as a convolution of a ^-correlated random process a{x), 
with {a{x)) — 0, and the modulation function G{x), 
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FIG. 2. (Color online) The inverse backscattering lengths 
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(6), (AS) 



(gray curve) and 1/L 



(6), (SGS) 



(black curve) for 



the second mode, n = 2 are shown; see Eqs. m and ([6|. The 
corresponding power spectra are calculated directly from the 
rough boundaries shown in Fig. [l] Regions of dominating 
SGS and AS mechanisms are shaded. 



From Eqs. ( |To| ) and ( [TT| ), one can readily recognize that 
the algorithm for constru cting a rough surface S,{x) with 
power spectrum ^ reads l^^*^ 



°° dx' , ,. sin(2fc+x') - sin(2A:_a;') 
a(x — X ) 



27r 



(fc+ - fc_)i/2a;' 



(12) 

Instead of a continuous a{x) we used a discrete series 
of random numbers to generate 60 discrete ^(x)-values 
at a; = mL/60 {m = 1, 2, 3, . . . , 60). The continuous 
roughness ^(x) was recovered by spline interpolation. 
After such a procedure the resulting boundary profile 
is no longer a pure gaussian random process and its 
power spectrum W{kx) is not perfectly rectangular as in 
Eq. ([9]). However, it still approximates the peaked struc- 
ture with the same fc^^-range. Note that the boundary 
has mainly two scales: 27r/(/c+ + fc_) (fast oscillations) 
and 27r/(fc+ — fc_) (beating oscillation) which can both 
be seen in Fig. [T] 

The inverse backscattering lengths, l/Ln'^'^'^^^ and 
1/Ln^''"^'^^\ following from SST for the constructed 
boundary profile are shown in Fig. [2] for the second mode 
n=2 (similar behavior for the first mode n = 1 is not 
shown) . We see that the inverse AS backscattering length 
([5]) shows the expected peak at fc„ « 75 m~^, determined 
by W(2kn), whereas the corresponding curve for the in- 
verse SGS backscattering length ^ dominates at much 
lower frequencies due to the influence of S{2kn). Note 
that as kn — > both curves in Fig. [2] diverge due to the 
factor k~'^. However, these divergencies are characteris- 
tically different for small frequencies. For the AS mech- 
anism the factor VF(2fc„) vanishes, limfe^_j.o W{2kn) = 0, 
while for the SGS mechanism the factor S'(2fc„) remains 
finite, limfe„_^o S'(2fc„) = 87r(fc5 - k^_)/5{k+ - k_f > 0. 
We conclude that the gap in the waveguide transmission 
for small values of fc„ is thus primarily induced by the 
SGS mechanism, unlike the transmission gap at around 
kn — 75 m~^ which is associated with the AS mechanism. 

In Ref. 37 we have demonstrated that the type of 



experimental setup we use (see Fig. [I]) allows for mea- 
surements of the scattering matrix Snm (in the mode- 
representation). Here we focus on the squared modulus 
of the diagonal elements of the scattering matrix which 
are proportional to the transmittance (jsjlj^^iSL 



{Tn 



(13) 



The properties of the antennas and their coupling to the 
waveguide modes are captured in the coefficients c„. In 
the experiment they are not directly accessible, but can 
be determined by fitting the transmittance through the 
corresponding mode for a roughness-free waveguide. A 
theoretical study of antenna coupling^^ predicts that the 
Cn depend only weakly on frequency such that we can 
regard them as fitting constants. Due to an imperfect 
coupling to the antennas and losses at open ends, the 
maximal measurable transmission is of the order of 10~^ 
as reflected in the values of ci — 6.33 x 10"^ and C2 — 
3.33 X 10^ (for details see Ref.f57|). 

The resulting transmission curves as measured for the 
first and second mode feature two transmission gaps [see 
gray curves in Figs.jSja) and[3]^b), respectively]. As pre- 
dicted by SST (see solid black curves for comparison) , the 
first transmission gap is located near the opening of the 
modes and the second gap appears around kn = 80 m~^. 
This clear correspondence between the experiment and 
theory allows us to associate the observed gaps with the 
SGS and AS mechanisms, respectively. To further prove 
this association we also calculate theoretical curves where 
the SGS term is explicitly neglected [see dashed curves 
in Figs. [3]^ a) and[3]^b)]. These curves only contain the 
AS terms and, indeed, do not display the previously 
identified SGS-gaps. We thus arrive at the central mes- 
sage of our Rapid Communication: transmission through 
surface-disordered waveguides manifests clear signatures 
of the AS and the SGS mechanism. The SGS contribu- 
tion to transmission may even dominate the conventional 
AS contribution and both mechanisms can introduce ex- 
tended band gaps in the transmission. 

SST thus reproduces all experimental features qualita- 
tively correctly. There are, however, systematic quantita- 
tive discrepancies between experiment and theory. Most 
prominently, the frequency interval of the observed trans- 
mission gaps are shifted to higher frequencies than those 
predicted [see Fig. 3]Ja) and[3f^b)]. In particular for the 
second mode this shift is quite sizable. To check the ori- 
gin of these discrepancies we performed detailed numeri- 
cal simulations for the same scattering setup. We calcu- 
lated all scattering matrix elements and the correspond- 
ing scattering wave functions (see Fig. |4]) in the waveg- 
uide with an extended version of the modular Green's- 
function methocpSIMI which is based on a finite-difference 
approximation of the two-dimensional scattering geom- 
etry. Both the AS gap and the SGS gap for the nu- 
merical transmittance T2 match the experimental values 
[see dashed gray curve in Fig. [sjc)]. We thus conclude 
that the observed shifts and broadenings of the gaps are 
not experimental artifacts, but rather details omitted by 
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FIG. 3. (Color online) Measured transmittance T„ (gray 
curve), for (a) the first (n = 1) and (b) the second mode 
(n = 2). The corresponding theoretical prediction is shown 
as solid black curve [see Eq. ([3|]. For the dotted line the SGS 
is neglected explicitly. The shaded areas indicate the gaps 
caused by the dominating scattering mechanism, (c) Numeri- 
cally calculated transmittance for the second mode T2 (dashed 
gray line), compared to measured transmittance (solid gray). 
The markers indicate the positions of the numerically calcu- 
lated wave functions shown in Fig. |4] (d) Numerically cal- 
culated transmittance (dashed gray curve) compared to the 
theoretical prediction (black solid curve) for reduced disorder 
strength a (50% of experimental value). 



SST. Since the latter was derived for perturbatively small 
disorder amplitude a, an obvious source of the theoret- 
ical shortcomings may be the sizable value of a in our 
experiment which goes beyond the first-order correction 
inherent in the AS mechanism as well as the second-order 
correction in the SGS mechanism. To test this assump- 
tion we compare the numerical calculations and the corre- 
sponding theoretical predictions for a waveguide system 
with the same boundary profile, but with half the rough- 



ness value of a [see Fig. [Sj^d)]. In this case we, indeed, 
obtain much better agreement between theory and nu- 
merics: The AS gap is now both less broad and shifted 
to the predicted interval (note that gap positions do not 
depend on a in SST) and the SGS gap is suppressed, as 
its dependence ^ cr^ in Eq. ^ would suggest. 
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FIG. 4. (Color online) Numerically calculated wave functions 
for transmission in the second mode. The corresponding k2- 



values of (a) 23.46 m 
indicated in Fig.jsJ^c). 



(b) 64.02m~\ and (c) 84.39m"^ are 



Our numerics also allow us to visualize the scatter- 
ing wave functions at characteristic frequency values (see 
Fig. |4|: at fc2 = 23.46m~^ the second mode is fully 
backscattered due to the SGS (a). Once this SGS gap 
is passed full transmission is recovered (b), followed by 
the AS-gap (c) where complete backscattering is found. 

To summarize, we investigated coherent transmission 
through a quasi-one-dimensional waveguide with surface 
disorder. Using a microwave measurement setup we 
find pronounced transmission gaps in predetermined fre- 
quency intervals, each of which we can associate with 
a specific surface scattering mechanism, i.e., the ampli- 
tude scattering and the square-gradient scattering. An 
observed shift of the amplitude scattering gap could be 
attributed to the non-vanishing disorder strength. Our 
investigations show that even relatively short waveguides 
can exhibit effects predicted for systems with long-range 
correlations leading to drastic changes in their transmis- 
sion properties. This result opens up the way to design 
boundary profiles to induce desired transmission proper- 
ties in surface-disordered scattering systems. 

Support by the DFG within the research group 760 
"Scattering Systems with Complex Dynamics" is ac- 
knowledged by O.D., H.-J.S., and U.K. F.M.L and 
N.M.M. gratefully acknowledge the support of the SEP- 
CONACyT (Mexico) Grants No. 80715 and No. 166382. 
J.D., F.L., and S.R. acknowledge support by the Vi- 
enna Science and Technology Fund (WWTF) through 
Project No. MAG9-030, by the Austrian Science Fund 
(FWF) through Project No. P14 in the SFB IR-ON, and 
computational resources by the Vienna Scientific Cluster 
(VSC). 



* lotto. dietz@physik.hu-berlin.de| 



lulrich.kuhl@unice.frl 



5 



^ A. A. Maradudin, Light scattering and nanoscale sur- 
1^ face roughness (Springer, New York, 2007) ISBN 
9780387356594 0387356592. 

^ F. G. Bass and I. M. Fuks, Wave scattering from statisti- 
cally rough surfaces (Pergamon Press, Oxford; New York, 
1979) ISBN 978-0080198965. 

^ G. V. Frisk, Ocean and Seabed Acoustics: A Theory of 
Wave Propagation (Pearson Education, Prentice Hall, Up- 
per Saddle River, NJ„ 1994) ISBN 9780132441544. 
E. I. Chaikina, S. Stepanov, A. G. Navarrete, E. R. 
Mendez, and T. A. Leskova, Phys. Rev. B, 71, 085419] 
[~ (2005) 

° M.-C. Phan-Huy, J.-M. Moison, J. A. Levenson, 
S. Richard, G. Melin, M. Douay, and Y. Quiquempois , 



Journ. of Lightwave Techn., 27, 1597 (2009) 



K. K. Lee, D. R. Lim, H.-G. Luan, A. Agarwal, J. Foresi, 
and L. G. Kimerling, Appl. Phys. Lett., 77, 1617 (2000) 
A. J. Maker and A. M. Armani, Opt. Lett., 36, 3729 
(2011). 

A. Chutinan and S. Noda, |Appl. Phys. Lett., 85, 3739| 



(19 99) 

■^"TT Roberts, F. Couny, H. Sabert, B. Mangan, T. Birks, 



J. Knight, and P. Russell, Optics Express, 13, 7779 (2005) 
^° L. Gan and Z. Li, Front. Optoolcctron., 5, 21 (2012)^ 
A. A. Maradudin, ed., Structured Surfaces as Optical Meta- 
materials (Gambridge University Press, Gambridge, NY, 
2011) ISBN 978-0-521-11961-a 
G. Fishman and D. Galeck, 
(1989)1 



Phys. Rev. Lett., 62, 1302 



^A. E. Meyerovich and A. Stepaniants, [Phys. Rev. B, 60, 
r 9129 (1999) 



and 



S. Zhan g, P. M. Levy, and A. Fert. jPhys. Rev. B, 45, 8689 
I (199 2)1 

F. Gamiz, J. B. Roldan, J. A. Lopez- ViUanueva, 
P. Cartujo-Cassinello, and J. E. Garceller, ' Appl. Phys.| 
I Lett., 86, 6854 (1993) 

^T. E. Huber, A. Nikolaeva, D. Gitsu, L. Konopko, 
M. J. Graf, J. Appl. Phys., 104, 123704 (2008) 

" G. B. Akguc and J. Gong, Phys. Rev. B, 78, 115317 (2008). 
J. Feist, A. Backer, R. Ketzmerick, S. Rotter, B. Huck- 
estein, and J. Burgdorfer, | Phys. Rev. Lett., 97, 116804 

|_ (2006)1 J. Feist, A. Backer, R. Ketzmerick, J. Burgdorfer, 
and S. Rotter, Phys. Rev. B, 80, 245322 (2009) 

^® G. Maire, F. Drsek, J. Girard, H. Giovannini, A. Talneau, 
D. Konan, K. Belkebir, P. C. Chaumet, and A. Sentenac, 
[Phys. Rev. Lett., 102, 213905 (20"09)| 



^° K. Ferry, S. M. Goodnick, and J. P. Bird, Transport 
in Nanostructures, 2nd ed. (Cambridge University Press, 
Cambridge, UK, 2009). 

J. B. Khurgin, Y. Dikmehk, P. Q. Liu, A. J. Hoffman, 

M. D. Escarra, K. J. Franz, and G. F. Gmachl, AppL] 
' Ph^s. Lett., 94, 091101 (2009) 

M. Y. Han , B^O zyilmaz, Y. Zhang, and P. Kim, Phys] 
I Rev. Lett.798,~206 805 (2007) 
" E. R . Mucciolo, K. H. Castro Neto, and C. H. Lewenkopf, 

Phys. Rev. B, 79, 075407 (2009) 
'^'^ M. Evaldsson, I. V. Zozoulenko, H. Xu, and T. Heinzel, 

Phys. Rev. B, 78, 161407 (2008) 
25 J. Feilhauer and M. Mosko, arXiv:1203.6512v2 (2012) 
'^^ V. V. Nesvizhevsky, H. G. Borner, A. K. Petukhov, 

H. Abele, S. Baefiler, F. J. Ruefi, T. Stoferle, A. West- 

phal, A. M. Gagarski, G. A. Petrov, and A. V. Strelkov, 

Nature, 415, 297 (2002) 

T. Jenke, P. Geltenbort, H. Lemmel, and H. Abele, [Nat] 
Phys., 7, 468 (2011) 

A. I. Hochbaum, R. Chen, R. D. Delgado, W. Liang, E. C. 
Garnett, M. Najarian, A. Majumdar, and P. Yang, Na- 
ture, 451, 163 (2008). 
2^ W. Zhao and J. W. Ding, Europhys. Lett., 89, 57005 
(2010). 

^° F. M. Izrailev, N. M. Makarov, and M. Rendon, Phys. 
Rev. B, 72, 041403(R) (2005). 

F. M. Izrailev, N. M. Makarov, and M. Rendon, Phys. 
Stat. Sol (b), 242, 1224 (2005). 
^2 M. Rendon, N. M. Makarov, and F. M. Izrailev, Phys. 
Rev. E, 83, 051124 (2011). 

M. Rendon, F. M. Izrailev, and N. M. Makarov, Phys. 
Rev. E, 84, 051131 (2011). 

F. M. Izrailev, A. A. K rokhin, and N. M. Makarov, PhysT] 
Rep., 5 12, 125 {2QVT^ 



34 



MrRendonTP^Mrizrailev, and N. M. Makarov, Phys. 
Rev. B, 75, 205404 (2007). 
'■^^ U. Kuhl, F. M. Izrailev, and A. A. Krokhin, Phys. Rev. 
Lett., 100, 126402 (2008). 

O. Dietz, U. Kuhl, H.-J. Stockmann, N. M. Makarov, and 

F. M. Izrailev, Phys. Rev. B, 83, 1342 03 (2011) 

T. Tudorovski y, R^ Hohmann, U. Kuhl, and H.-J. 

Stockmann, J. PhysrA, 41, 275101 (2008) 

S. Rotter, J.-Z. Tang, L. Wirtz, J. Trost, and 

J. Burgdorfer, Phys. Rev. B, 62, 1950 (2000). 

F. Libisch, S. Rotter, and J. Burgdorfer, "Edge scattering 

in graphene nanoribbons," arXiv:1104.5260[New J. Phys. 

(to be published)]. 



